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ABSTRACT 


A study of single-span initially bent columns is made and extended 
to include two-span columns. The axial force in these columns results 
from heating (or an equivalent) and from bowing. A snap-through 
phenomenon is considered in the analysis of the single-span and 
observations are made concerning snap-through in the two-span column. 
A digital computer program is developed to carry out the analysis and 
a description of the program enables its use by anyone familiar with 


the FORTRAN language. 
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ieee T RODUCTION 


A. GENERAL REMARKS 

The subject of initially bent columns and shallow arches is 
treated extensively in the literature. Hsu [Ref. 6] studied the 
dynamic stability of shallow arches against snap-through; Hoff and 
Bruce [Ref. 5] analyzed dynamically the snap-through of shallow 
arches; Fung and Kaplan [Ref. 2] studied the buckling of low arches 
or curved beams of small curvature; Timoshenko [Ref. 8] has also 
discussed the problem; Huddleston [Ref. 7] studied the effect of 
axial strain in elastic columns; Gjelsvik and Bodner [Ref. 3] 
approached the problem on the basis of erent criteria. There 
have been other similar studies made. However, these treatments 
of the problem have all involved only a single span. Brock [Ref. 1] 
studied the single span using Fourier analysis in such a way as to 
provide for assembly into a multispan continuous column, and gave 
conditions for the solution of multispan cases. However, he did 
not actually arrive at numerical solution in any particular case, 
only suggesting how the computations might be carried out. There- 
fore, it was the objective of this thesis not only cOustudy Lie 
single-span column in greater depth, but to develop its analysis 
in a manner such that it could most easily be extended to a multi- 
Span problem, and actually to make this extension to the case of a 


two-span column. 





Bey OCOPE OhaWORK 

The theory of the single-span Boil etl is first developed for a 
column which is initially nearly Baan. This development takes 
into consideration an axial strain induced in the column and also 
allows for the external application of end moments. It is assumed 
that the configuration, with and without load, lies in a single 
plane. The theory is then extended to include the two-span column. 
A digital computer program is developed to carry out the analysis, 
and a description of the program enables its use by anyone who is 
familiar with the FORTRAN language. 

Analysis beyond the two-span problem is not carried out in this 
thesis, although it 1s envisioned that multi-span problems can be 


investigated in much the same manner as described herein. 


ee O LAT LON 


A - cross sectional area 

a - Fourier coefficients for initial configuration 
b - Fourier coefficients for final configuration 

E - Young's modulus of elasticity 

I - moment of inertia of column 


Peo elenetii ol Span 


M - left end moment 


N - right end moment 

n - term index in Fourier analysis 
P - axial compressive load 

eg Lider asenritical buckling load 
Q - left end reaction force 


R - right end reaction force 





lateral load in two-span column 
distance from left end of span 
imieieaeeder lection 
im Taaerleat.on 
g 
ieee train 
thermal unat strain 


curvature 


excess of arc length over chord length in initial 
configuration 


excess of arc length over chord length in final 
configuration 


6) le) 7a 


Cideslepe OL Span 





Meee levee bOrMeENT OF THEORY 


A. SINGLE SPAN 


The study of the single span involved an initially bent column 


loaded as shown below: 





Figure 1. Loading on single span column 


The end supports are frictionless hinges which do not move from 
their original positions. The moment loading shown is applied by 
external agencies. The initial, no-load configuration of the column 


is presumed to be given by the equation 


= . Nt xX 
ve =o. aa Sin Gy A (1) 


where here and later the summation is with respect to n and goes from 
n=l to n=. For computational purposes the upper limit is taken to 
be a finite number, which will be discussed in a later chapter. 

Next suppose that the column is uniformly heated so as to provide 
thermal strain. (It should be noted that an equivalent situation 
would result from moving the ends of the column toward each other. 
However, throughout, we will speak in terms of heating the column 
rather than displacing its ends.) The — configuration can be 


presumed to be represented by the equation 


=f nin oeees (2) 
y Do sin (2), | 
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The immediate purpose is to find the b's and thereby the new 
configuration. 
There is an axial compressive force P set up in the column such 


nets tne wntiet Strain is 


€ = et = eae C3) 


where “or is the thermal strain, A is the cross section area, and E 
is Young's modulus of elasticity. (Throughout, we presume that Hooke's 
law is satisfied.) 

The excess of arc length over chord length can be adequately 


represented for the initial oe by the approximate relation 


2 
Tl 
ay * dx = Ena? (4) 


Likewise for the new Sn ee 


12 


Dect oy 20 5 
age 7 (5) 


Thus, the change in length due to heating is 


2 
oy En2(b * - a 2) = z 6 
5 . (7) n (b_ a ) Le,. PL/AE (6) 


The bending behavior of the column must also be accounted for. By 
static equilibrium 
Q = -R = (N-M)/L. (7) 


The bending moment at any point is given by 
B.M. = QX + M - Py, (8) 


and this accounts for the change in curvature 


ay oe : | 

™ 2 2 - {narx 

A —— | 7 == _—)} ok ae 

ao Ge 2 (* n (b_ a.) sin (2%) ’ (9) 





According to the elementary formula, 
io Uign =: 1th aes (10) 


Thus, with a little manipulation, we find 


>: 2 2 _ : igs 2 
[n P Ob, a_) Pb sift ( ae Oe 4 =0. (11) 


where we have introduced the notation 
Po insti nee (12) 


This is the Euler critical load which appears naturally in this 
Se lemiattone 10 deal successfully with Bq. (11) it is useful to 


express Qx and M as Fourier series in the range 0 < x < L. 


Cement) Sy" Gyn) sin (nx (13) 
= (2M/nm) Yf1 - (-1)") (1/n) sin (==). (14) 


Thus, 


7.23 cyt (208), pee) "(2 ‘ye n@P (b_-a.)-Pb_ sin (22), (15) 


which will be satisfied if each of the terms in braces vanishes. 


Accordingly, we set 


_)}_2 men) 200 7 ee n= 11,2508. (16) 
Dy Sis 22 i») { 2 }+ - =): om (2 ale rp B) i 


However this equation still involves the unknown, P. Equation (6) can 


bewwrertten in the form 


P = Ake_ -n2AE/H4L2 £ n@(b 2- a 2). Gigs 
T n n 


Rewriting Eq. (17) once more in the form 
A ce Oe (18) 


where 


10 








£(P) = PL/AE + (12/4L)En?2 (b_? E a *) mie (19) 


T? 
we seek the zero (or zeros) of the Periacee function f(P). Because 
of nonlinearity, computational difficulties were encountered in find- 
ing the significant root. Of several methods attempted, the most 
Suitable was the Newton-Raphson iteration scheme. A study of the 
curve f(P) revealed that the method would usually converge on the 
first root if the iteration was started at a value just slightly 
emailer than Pa 

Hmaemvaite or P thus being known, Eqs. (16) can then be solved for 
ime coefficients bo 

Eieimee2 shows the initial and final deflection curves of a column 
computed using this analysis. The final deflection shows the column 
after having undergone thermal strain. The end moments were zero for 
this simple case. 

Further analysis of the single span problem will be given in 


Section III of the thesis. 


be WO SPAN 

We now consider two such spans, the column being joined at the 
intermediate support so as to maintain continuity of slope and moment 
Hee iimdamoredetlection. In general, ieee: axial loads wiil be 
acting in the two spans. As before, the problem is to determine the 
final configuration after heating, the initial configuration having 
been specified. 

The obvious method of analysis is, in effect, to consider two 
adjacent single spans, each of which is as shown in Fig. 1 and is 
analyzed as in the preceeding subsection. The adjoined spans are as 


shown in Fig. 3. Separating them and considering their mutual 


ial 





NIVYLS IWWYSHL ONIOOYSGNN NNNIOD NYdS@STIONIS °2 SYNSIS 


Co 
SSHON! - HLONST NvINI09 
08 | ool Ov! 02! 00! 08 09 Ov 02 

: : +0 
i O 
mM 
ma 
_ 
m™m 
> 
sa 
NOILOSTS3G WILIN! 22 

t 


iam, S 


mA v 
NO1103914530 WNI4 


SSHON | 





interaction and the continuity at their junction, we have (see 


Fig. 4) 
ine te (20) 
6 = 6 
Pi 2 cy) 
If the common value Ny = M can be determined, the problem is 


solved in the sense that each span can be now analyzed separately. As 


Baeilatter of fact, in the solution for Ny = M, > we Will actually have 


determined Py and P. and the Fourier coefficients Dis for each span. 


M, Nz 


Figure 3. Loading on two-span column 





Figure 4. Separation into two single-span columns 


The left and right end slopes of each span may be found by 


evaluating dy/dx at x = 0 and x = L, respectively. 








nb 
@ = b eb. (22) 
n tb 
os aD) ; ae (23) 
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A method of solution for the two span problem is outlined as 
follows. First assume a common value for Ny and M,: (This will be 
called the junction moment.) Then analyze each span separately, using 
the analysis of the preceding section, and determine the slopes at 
the junction. Compare the two slopes and find the difference between 
them, which will be called error. Next repeat the above for another 
value of moment and find the error. A plot of error versus moment 
Will indicate the value of junction moment that will cause the error 
(in slope) to be zero. At this point, correct values for the Fourier 
coefficients be will be known for each span so that the final config- 
uration can then be computed and plotted if required. 

If the plot of slope error versus junction moment is nearly 
linear, a slope intercept method of iteration may be used to determine 
where the error goes to zero. This method simply takes two points, 
gets the equation of the line passing through these points, and then 
solves for the zero of this linear approximation.. This 1s repeated, 
uSing a new point and one of the old points, until the desired 
convergence is reached. If the plot of slope error versus junction 
moment is significantly nonlinear (this occurs if the axial compressive 
load is a large fraction of the Euler load), then the slope intercept 
method is no longer satisfactory and a method such as successive 
decimation must be used. In this context, decimation means recalcu- 
lating over an interval of interest, using increments one tenth as 
large as previously. The program decimates four times so that the 
final increment is .0001 times the initial increment. This method is 
well known and will not be described here. The accuracy of results 
obtained by using this method depends upon the size of the increment 


used in stepping. The smaller the increment, the greater the accuracy. 


14 





Sitewgucstiom bere is where to start the iteration. If a good 
approximate figure is known, then the increments may be small. 
If not, the increments must be larger. Unfortunately there is no 


quick way of determining a good starting approximation. 


eee INITIAL DEFLECTION 

For the single span problem, any initial configuration may be 
specified subject only to the condition y(0)=y(L)=0. If the 
configuration is specified in the form of a polynomial in x, the 
axial distance from the left end, the procedure to be described 
later quickly determines as many of the Fourier coefficients a, as 
may be desired. 

However, for the two span problem the necessity for continuity 
places extra requirements on the initial deflection. Not only must 
the deflection be zero at each end of the spans, but the slope and 
curvature must be eontinuous at the junction. Accordingly, such a 
two span configuration is obtained as follows. An arbitrary config- 
uration (conveniently in the form of a polynomial) is assumed which 
has zero displacement at the left of span one and at the right of span 
two. (see Fig. 5). Then, regarding this as a simply supported beam, 
sufficient lateral load is applied at the ansbeiemed lense support to 
cause the displacement to vanish at this point. The corresponding 
deflection curve is calculated and its superposition with the 
originally assumed configuration gives a result which now has zero 
deflection at each of the three support points. In general, the 
column will thus start out with pa ree moment except at 


the extreme ends. 


Age 








Figure 5. Original specified configuration 


Pigure 5 illustrates the originally specified configuration 
before the intermediate constraint is enforced. It is not dif- 
ficult to show that the necessary lateral load is given by 
3(L,+L,) EIA 
Force = W = ——y—~—y-—-—~_—_ downward. (24) 


os 


It can also easily be shown that the deflection of a beam due to a 
downward force W at point B is represented by 


WL, . ; : : Cree 
: Ao ao 2 ee 25 
6E1(L,4L,) (Lith, )° x-x - Lo°x + i | (25) 


where the < > symbol has the meaning 


0 Tae aes <4, 
Soe 2 i = 


(x-L,)° if x OL 


Thus, by superposition, the new configuration is 


WL, 5 : . (L,+L,) on 
E : Mee eS cy 2 Lo C26 
anne x) CEI, 4h) (L,+L,) Oe ne 1 


For the first span the expression for the configuration is 


WL, ps4 : : ‘ Can 
= - - - L " 
yee) GEI(LGEL) (L,+L,) 2 0 
jee 
where 0 < x < L,- 
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For the second span 


WL. Cx*4+L_ ) ieee) 


a fe 2 a 2 of i 2 i 2 “s 5 
ca fx +L, ) SET(Lo aE) (L+L,) (x +L, ) Ly + T(x y (28) 
ie 2 i 
where x* = x-Ly is the distance along span two from its own left end. 


Now the Fourier coefficients can be found for each span. If f(x) 
is a polynomial, the program described in Appendix A makes the deter- 
Mination. Thus we have a deflection satisfying the necessary initial 


conditions for the two span problem. 
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ie Discuss lon 


ie ANALYSIS OF THE PROBLEM 

In analysis of the single-span problem, the effects of the 
application of end moments were investigated. This means that 
end moments were externally applied at one or both ends of the 
column with or without the presence of axial strain. It was 
found that a snap-through phenomenon could be observed using 
this analysis. That is, by increasing one or both end moments, 
the configuration could be seen to "jump" from one configuration 
to a distinctly different configuration with only infinitesimal 
variation in moment. Figure 6 illustrates this phenomenon by 
showing a column that is first heated up, and then subject to a 
left end moment. Curve 1 represents the initial configuration 
of the column. Curve 2 is the configuration of the column after 
undergoing thermal strain. The remaining curves show the config- 
uration under the action of several different values of end moment. 
It is seen that the first few 1000 inch pound increments of applied 
moment cause only slight changes in configuration, but that the 
change from 6000 to 7000 causes a relatively leas change which 
has the nature of a jump, or snap-through phenomenon, following 
which addition of increments of applied moment again result in 
only minor adjustments of the configuration. 

Figure 7 further illustrates the phenomenon by showing a plot 
of centerline deflection versus aoe where time is used in the 


sense that first thermal strain is inereased uniformly, with no 


18 





NWN109 NVdS-SIONIS NI NONAWONSHd HONOYHI-dVNS °9 3YNDIS 





°a1 “NI 0006 
Sane -Nieoooe 
"G1 "NI 0002 


SSHONI - HLONAT NWA100 


O91 00! 08 09 





SSHON | 





NO!LVYEND 1 SNOO 
WILINI 


"a1 "NI 0009 = 

°a1 °NI 000S = 

"a1 °NI O00r = 

O=W 

INVLISNOO = NI/NI 100°0 ="93 


men @ 1 Os 1340 


Lo 





CENTRE eS NOU ESIEIGE 2G) ro BUNS SC) sails Sent hie Se USNS 6 seh s 


BONO L=aVNS daly 
SNOILVYEIA MOHS 
GINOM S3SVO WALOV 





ANAWON 
LNSWOW GN3 ONISYSHONI ‘SYNLVYSDWSL LNYVLSNOO ONS ON ‘ONILW3H 
ee ) LNAWON i tp 


LNANSYONT TIVWS 

JNO AG LN3WOW 40 

i Sov saGeN | wok ad aov Id 
SSxNVL HONOYHL-dYNS 


LNAWSYON | 
LNAWON 
INO 


NO!103174350G YsALN3a0 


20 








end moment, and then thermal strain is held constant and the end 
moment iS increased at a uniform rate. 

It may be noted that a similar analysis could be made for a 
column under lateral loading instead of end moment loading. 

It will be of interest to explain why the phenomenon of snap- 
through is considered to be of great importance. In a piping system 
which is heated up, considerable energy is stored in the form of 
Strain energy. If snap-through or a similar phenomenon could take 
place, some of this energy might be suddenly released and have a 
catastrophic effect on the piping system or adjacent installations. 
In the case of a single-span straight-line system, if the ends 
should be pinned, so that no moments may be applied, it is clear 
that snap-through cannot take place. If the ends are fixed to 
heavy equipment, such as pumps or vessels, the restraining influence 
of these items would indeed result in end moments being applied to 
the pipe, but it is difficult to envision these moments actually 
resulting in snap-through. 

However, in the case of multi-span straight-line configurations, 
at the intermediate supports adjacent spans exert moments on each 
other, and it is not immediately clear that as the system is heated, 
perhaps nonuniformly, these intermediate ae Mihi not amercase 
in such a way as to cause snap-through. The intuitive feeling that 
this would not take place was insufficient basis for asserting that 
it could not happen. Accordingly, considerable attention was focused 
on this possibility. It was only after studying many different ube atk eulra 
configurations and thermal loading sequences shee a reasonable under- 


standing was achieved of the restraining influence which one span has 


PAL 





upon an adjacent span and of how this serves to prevent the phenomenon 
of snap~-through, even if one end of a waisdl ROSTER configuration is 
subjected to increasing external moment. 

alee rilyy. the simplest multi-span configuration is that having 
two-spans. Reson anpigs this case was studied in detail and computer 
programs were written for its analysis. One of the cases which was 
studied is illustrated in Fig. 8. In this case a negative moment 
of increasing magnitude was applied at the extreme left causing the 
column to go through the succession of configurations shown. (Inci- 
dentally, although initially there was axial compression in the column 
due to initial heating, by the time the last configuration shown was 
obtained, the axial load in the left span had become tensile.) 

This example provided the insight to see why snap-through, at 
least in any sense similar to that for the single-span column, is 
not possible for two-span columns. The original, incorrect, view 
was that moment esate at the left end of a two-span column was, 
in effect, simply a way of applying moment to the left end of the 
right-hand span; increasing the leftmost moment sufficiently would 
eventually apply a very large moment, at the junction point, to the 
Hetieend of the right span. This) 1s true enough, but the continuity 
between left and right spans does not permit the latter to snap- 
through as would be the case if it were the only span and the 
junction moment were directly applied at its left end. In effect, 
the left span acts as a stabilizing influence to prevent snap-through 
Ciena right Span. 

The situation may be visualized more clearly by imagining that 


the continuity between the two spans is "one-way" in the sense ona 


ZZ 
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the mechanism shown in Fig. 9. If the two spans were joined in this 
manner, snap~through of the right span would indeed take place when 
the junction moment reached a critical value. However, with the 


true continuity of a continuous column such behavior can not take 


race . 
WITH THIS PIN REMOVED, AS 
SHOWN, CONTINUITY IS 

M ONE NVAY 





Figure 9. "One-way" continuity of two-span junction 


Attempts to create a situation in which snap-through takes place 
in a two-span column have proved unsuccessful. Study of the behavior 
in these cases has added strength to the conviction that such a 
phenomenon cannot take place, even under the action of increasing 
externally applied moments at the ends or at the intermediate supports. 
However, this growing conviction cannot be understood as a convincing 
demonstration that the phenomenon cannot indeed take place, perhaps 
only for certain unfortunate combinations of physical parameters. 
Accordingly, further study of snap-through in multi-span columns is 


suggested as an appropriate topic for further investigation. 
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B. FOURIER REPRESENTATION OF THE CONFIGURATIONS 

A question arises of the ability of Fourier sine series, as 
employed herein, to represent non-zero curvature at the ends of 
Spans where each sine term vanishes. Certainly, the value of 
curvature given by such a series must be zero when x=0 and x=L. 
If we have initial deflections which have non-zero curvature at 
the ends of the span (or spans), our Fourier series must be able 
to represent this adequately. ep if we have non-zero moment 
at these points (due to the application of external moments or due 
to enforcing an intermediate constraint) this corresponds to non- 
zero curvature change, which must be adequately represented by our 
Fourier series. 

Accordingly, actual plots of curvature, as determined from 
the Fourier series, were made in several cases where the value of 
curvature could be determined by other means. Figures 10, 11, and 
12 show a comparison of curvature in a case where the initial 
configuration (in the sense of Fig. 5) had uniform negative 
curvature and where the intermediate constraint produced a piece- 
wise linear moment change (i.e., curvature) that was zero at each 
end and a maximum at the constraint. The final curvature is the 
superposition of these, consisting of two straight segments rising 
from -.001332 to .000666 and back down to -.001332. (In Figs. 10, 
11, and 12 these straight line segments show a very small waviness 
due to the method of plotting.) These figures also show the Fourier 
series representation of ee using 25, SOmmondm OO0sterns , 
respectively. It is seen that although indeed the Fourier repre- 


sentation of curvature does vanish at the constraint points, this 


25 
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is a very local phenomenon and elsewhere the Fourier series provides 
a good representation of the correct values. 

Another example is shown in Fig. 13 where values of curvature 
as represented by Fourier series are plotted for different values 
of left end moment in a case of snap-through. Fairing straight 
lines through the left ends of the curves gives intersections with 
the left vertical axis which correspond to the actual moments applied 
eaecthne left. 

Probably the most convincing verification of the adequacy of the 
Fourier representation is the relative insensitivity of significant 
results (in a particular problem) to the number of Fourier terms 
employed. In a particular case of snap-through, Table I below shows 
that snap-through occurred in the same moment interval for all values 
erm down to as Jiow as 3. It is interesting to note in Table II that, 
even for n=l, the axial load obtained is in very good agreement with 
the load obtained ee 100 terms. 


Table I. Values of end moment between which snap-through 
Gceurred, Lor different yalues of n. 


Moment (in.1lb. ) 


60450 


60460 


60470 





Note: Dots indicate the moment interval in which snap-through 
occurred. The tabulation is for a particular problem 
described on page 31. 
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tate aiuesOr axial load in pounds for different 
values of n and M. - 


Moment Cineiine ) 60450 60460 60470 


38489.000514 38491.906680 38489. 899526 
38489.000515 38491.906680 38489.899526 


38489 .000517 38491.906680 38489 .899527 











38489.011286 6679h, 0096523 38489 .906358 
38489.011298 38491.009624 38489 .906349 


The tabulation is for a particular problem described below. 





Due to computational difficulties in the vicinity of snap- 
through, values of moment taken at smaller than 10 inch-pound in- 


crements provided somewhat erratic results. 


The results tabulated in Tables I and II are for a steel pipe 
whose cross sectional area is 2.945 in*. The length is 200 inches, 
the moment of inertia of the cross section is 5.2 in’, and Young's 
modulus of elasticity is 30 x 10®psi. Moment is applied at only 
Suesenanor the span. The initial configuration was similar to that 


Sewn in fig. 6. 
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IV. ANALYSIS OF A TYPICAL TWO-SPAN PROBLEM 


Using the program described in Appendix A, configurations 
resulting from the application or many different combinations of 
end moments and strain on a two-span column can be obtained. The 
following cases illustrate several of the ways in which the program 
may be utilized. These examples involve a two-span steel pipe with 
cross sectional area of 2.945 in*. The left Span is 100 inches 
long and the right span 300 inches. The moment of inertia is 5.2 in! 
ema oung’s modulus of elasticity is 30 x 1a Psi. It should be 
noted, however, that a column of any constant cross section or 
material could have been used. 

Case I - As illustrated in Fig. 14, the pipe is subjected to 
thermal strain only. Curve 1 shows the initial configuration. 
Subsequent numbered curves represent the configuration correspond- 
ing to increasing values of thermal strain. 

Case II - Here the pipe is subjected only to an increasing 
moment. Figure 15 shows the sequence of resulting configurations. 

Case III - This example shows the resulting configuration of 
the pipe (see Fig. 16) with simultaneously increasing thermal strain 


and left end moment. 
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V. RECOMMENDATIONS 


If future studies are made uSing the analysis set forth in this 
thesis, the following areas are recommended. 

The problem should be extended to include multi-span columns. 
In the analysis of two or more spans it would be convenient to have 
a better way of obtaining the junction moment or moments. Perhaps 
a method can be developed that does not require an initial approxi- 
mation. 

Another area of study would be to develop a proof that snap- 
through cannot take place or develop a clear delineation of the 
circumstances under which it could take place. The study probably 


should focus initially on two-span columns. 
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Poet se DESCRIPTION OF COMPUTER PROGRAM 


A. GENERAL REMARKS 

The digital computer program used to carry out the analysis set 
forth in this thesis is a double precision FORTRAN language program. 
It 1s self contained except for one outside subroutine used to draw 
the configurations. This subroutine or a similar one should be 
available in any computer library. Various write statements are 
interspersed throughout the program to allow the user to follow the 
progress of the computations, and to help trace any errors. These 
statements may of course be left out of the program to lessen the 


amount of unneeded output. 


B. BASIC FUNCTIONS PERFORMED BY THE PROGRAM 

The basic program for the single-span plots the initial deflection, 
determines the axial load on the column, and plots the final deflection. 
The axial load is determined using a Newton Raphson technique on Eq. 
(19). Iteration continues until convergence is reached. 

For the two-span problem, a process 1s employed of repeatedly 
selecting values of the junction moment. Then .the two spans are 
analyzed and the two junction slopes are computed and compared. Then 
the correct value of the junction moment is found using the method of 
successive decimation. Although the slope intercept method is more 
accurate when it works, its unreliability causes it to be less 
desirable for use in the general program. After the junction moment 


is found, the final deflection is then plotted. 
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PLL IAG, DEFLECTION 


The Fourier coefficients for the initial deflection of the column 
are computed for the analysis previously described. The program 
accepts coefficients describing a polynomial curve which passes 
through the left support of span one and the right support of span 
two. Then the modification of this curve is computed which enforces 
the intermediate constraint; this modified curve is also in the form 
Pea polynomial, ~This yields two polynomial curves, one for each 
Mion, cr. Eqs. (27) and (28). Then the Fourier coefficients a are 
evaluated for each span. 

item oubroutine SINCOE 


This subroutine evaluates the integral 
1 


GA(m,n) -f a sin(ny)dy CA) 
@) 
mimes used in the evaluation of the Fourier coefficients. The 


method involves evaluating integral number 6, page 117, of Ref. Da 
See program listing in Appendix C. 
Peeivelvatlonwon Fourier Coefficients 
The Fourier coefficients are evaluated according to the 


following formula: 


L 


; | 
a= i fee sin (2) een 1 260. (A.2) 


0 
where L is the length of the span and n is the number of coefficients 


desired. For the left span, y, as given by Eq. (27), is substituted 
for the value of g(x), and for the right span Eq. (28) is used. 

The method carried out in the program can best be illustrated 
through the use of an example. Suppose the initial configuration f(x) 
is given by the equation 


f(x) = Ax* + Bx. (A.3) 
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Then Eqs. (27) and (28) become 


WL, 


= Ax2 a EE SESS ERO OR?) 
yes Axo + Bx SEY(L, tb) “ (L, +L,) x L, 


and 


WL, 


= oh Z oy i rhs 2 me eres 2 = 
y = ACx +L, ) + B(x +L, ) SET LiL, (x +L, ) [cyst (x +L, ) 


(L, +L) 3 | 
+ 


Lj (x*4+L, ) 
respectively. Equation (A.4) can be rewritten as 
y= Ax? + Bx - Cx + Dx3 


and Eq. (A.5) can be rewritten as 


= OG 2 ek = ze ae 3 Oh 
we A(x +L, ) + B(x +L, ) a6. +L, ) + D(x +L, ) Fx 


where 
WL(L, +2L,L,) 


Cc = 
3 
6EI(L,+L,) 


WL, 


D = SEL? 
6EI(L, +L, 


and 
_W 
Oi. : 


For the first span, Eq. (A.2) becomes 
Ly L, 


‘TT 
a = = a fx? sin (BE lax + a fs sin (27 Yen 
a 1 1 il 


0 
L, 


0 
Ly 
= (6 [> 2in( BP) ax + pf x3 sin (AE Je 
L L 
0 Ale 1 


0 





The integrals are all of the form 


L, 
fe cin (AEN ox 
a 
0 
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Cho 


(A. 


ce 


(A. 


(A. 


(A. 


uv) 


a S)) 


6) 


7) 


8) 


9) 


AG.) 


i) 


12) 





1 
Making the substitution y = — this integral may be put in the form 


1 
m+1 TT 
a m n+l i 
== y sin (ny Jay = 0 GA(m,n), (CAelion) 
0 
where 
es: 
ee. (A.14) 


and the definition, (A.1), is used. Thus Eq. (A.11) can be rewritten 


in the form 


AY) 
HN 
yee 


| %0n¢ 2,0 + BO*GA(1,n) - Ce“GA(1,n) + De%ea(34n)| 
1 


oa Aeon GAC i,n) CNS. 
Ly iL 


Note, from Eqs. (27) and (28), that the maximum degree of the 
adjusted polynomials is Max(3,m*), where m* is the degree of f(x). The 
number n of Fourier coefficients to be employed, and the integer m* 
are among the inputs to the program. Internally, m* is increased, if 
mecessary, to a Haan value of 3. Then SINCOE is used to generate 
EMmemmiumerical constants Gi(i,j), 1 = 1,2,06.., m*; j = 1,2, cane ie 
Thus, the coefficients are available for the evaluation of the Fourier 
coefficients as represented by such equations as (A.15), the coef- 
fmeted seneeb. 0, ..., in which are obtained in. an obvious manner. 

The Fourier coefficients for the second span are found in 
basically the same manner, except that now the computations are more 
numerous due to the necessity of expanding powers of (x*+L). It can 
be seen that for the second span Eq. (A.2) becomes quite lengthy since 


g(x) will contain many terms. Looking only at the first term in the 


example, we have, after expansion 
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p) 2 ly 
ase. 1 
an if A wee sin t — \dxoaoL x* sin ue dx 
2 L ih L 
0 : 0 
Le 
TT 
Fs fos (2E)o ee ; (A.16) 
al L 
4 Zz. 


irs can also be written in the form 


saecl ma 


Jl 


=) 2 
| > a olcncayn) + & | 


%ea(0n)]+ ie GAD 
a 2 


where now 


Ae) 
Il 
ca 
SIN 


D. VARIATIONS IN THE PROGRAM 

Different combinations of strain and end moments to be applied to 
the column can be incorporated into the program. This can be accom- 
plished easily through the use of the input cards described below. 
The program can also be altered through an additional DO-loop to show 
a "time history" of the column. The strain and/or end moments can be 
incremented in Rie! D0=i0Cp so that intermediate configurations may be 
plotted. This is useful in looking at snap-through. The program 


listed in Appendix C has this variation incorporated. 


Hee vor OF THE .PROGRAM 

This program may be used by inserting the following Al input ecards 
in the designated places at the beginning of the program. Each card 
contains only one entry. Note that these cards are made a part of the 
program deck and are not read in as data following the program. see 


program listing for example. All cards must be punched beginning in 


column seven. 


ud 





DIMENSION GA(MPOL,NFOR) - integer numbers must be 


MPOL 


NFOR 


NDEC 


Co = 


Co = 


CH = 


Cir— 


C2 = 


Ci= 


Ll = 


L2 = 


E 


W 


Z 


lejersdl 


lglensyZ 


SS = 


UU = 


XMOMT 


inserted for MPOL and NFOR 


= maximum degree of initial deflection function 


=eMmiUmMpomrormerolrier coetticients desired 


= number of decimations in determining junction 


moment 
iit dl fonction 
ihe hd Uunei Lon 
iiietal. function 
initial function 
iitelaletunecLon 


inet fUunCLLon 


coefficient 


coefficient 


coefficient 


coefficient 


@eoprlcrent 


Cecnel ecient 


tense or Tertl span 


length of right span 


Young's modulus of elasticity 
cross section area of the column in inches * 


: es e e e 
Momemeeor Inept id Gmeene column in inches 


= applied strain in left span in inches/inch 


iT 


applied strain in right span in in/in 


applied left end moment in pound inches 


applied right end moment in pound inches 


= starting moment for successive decimation 


iTrerat on 


DELTA = initial increment for successive decimation 


Notes: 


iteration 


1. All inputs are floating point numbers (D format) except MPOL 


and NFOR, which are integer numbers. 


2. There must be an input card for each input listed, even if it 


1s zero. 


U2 





3. Each initial function coefficient must correspond to the 


power of x it is associated with. That is 


fe) = Cox” + Cox” 4 °CHx’ + C3x° + C2x* + Cix 
(A polynomial of degree greater than six can be accommodated only by 
making internal changes in the program: these should be rather 
obvious. ) 

4. The program is dimensioned for a maximum number of 50 Fourier 
coefficients and a maximum combined span length of 400 inches. (This 
limitation is related to the fact that deflections are computed at 
two-inch intervals and are stored prior to output.) If the user wishes 
to extend these limits, the dimension statements must be altered 
accordingly. 

5. Likewise the call statement for the DRAW subroutine may have 
to be altered to accommodate the particular problem in question. (For 
example, to provide a desired scaling.) 

6. Where to start the iteration for the method of successive 
decimation, i.e., selection of the input XMOMT, poses a problem. 
Experience and perhaps a trial run will best produce a suitable start- 


ie POlne. 
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APPENDIX B ~ FLOW CHART 
SHA 


DIMENSION 
DATA INPUT 


XCOUNT 


it 


Gay 


ACOUNT 
BCOUNT 


tou 
i ep) 
O © 


ACOUNT = ACOUNT + [|.0 





. 


XCOUNT = XCOUNT + 1.0, 


COMPUTE FOURIER 
COEFFICIENTS, a, 





/\ 7 B\ AS 


in 


SW 








AA 


COMPUTE LOAD 


eo PUTEMEOURTER COEFFICIENTS, b 
COMPUTE END SLOPES 





COMPUTE AND PLOT 
FINAL DEFLECTION 


COMPUTE END SLOPE ERROR 
SEV MEMHOO OR eoUCEES Sl VE 
DECIMATION 


BCOUNT = BCOUNT+ | .0 





ee) 


— 
> 
<u> 


=NDEC 
=2 


s) 
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APPENDIX C - PROGRAM LISTING 
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